Pion condensation in a dense neutrino gas 
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We argue that using an equilibrated gas of neutrinos it is possible to probe the phase diagram 
of QCD for finite isospin and small baryon chemical potentials. We discuss this region of the 
phase diagram in detail and demonstrate that for large enough neutrino densities a Bose-Einstein 
condensate of positively charged pions arises. Moreover, we show that for nonzero neutrino density 
the degeneracy in the lifetimes and masses of the charged pions is lifted. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) predicts that the 
properties of hadronic matter hugely depend on the 
baryon and isospin density. For example, it is found that 
for small baryon chemical potential such that there are no 
baryons around a Bose-Einstein condensate of charged pi- 
ons should arise if the absolute value of the isospin chem- 
ical potential (///) becomes larger than the vacuum pion 
mass (m n ) This has been shown using the chiral ef- 
fective Lagrangian [H, 0, 01 • Furthermore the existence of 
apion condensate was also demonstrated in_lattice QCD 
[2 H H 0, HI j with random matrix theory 
Nambu-Jona-Lasinio (NJL) mod el Jl Ol. \ 
with the 0(4) linear sigma model [IM,[l6[ and with string 
theory inspired holographic models of QCD [ijl [H, Gil ■ 

The question we would like to address in this article 
is under which conditions a pion condensate can be real- 
ized in a macroscopic system of particles. In order for a 
macroscopic system to be stable, it should be electrically 
neutral. Since the isospin chemical potential is equal to 
the charge chemical potential (/xj = /iq), requiring neu- 
trality puts restrictions to the value of the isospin chem- 
ical potential. 

At low baryon chemical potential (fis ^$ m p , here m p 
denotes the mass of the proton), the electric neutrality 
constraint forces up ~ so that pion condensation be- 
comes impossible [2(| (see also Refs. [U, HU). If how- 
ever a finite density of neutrinos which is in equilibrium 
with hadronic matter is present, these considerations will 
be modified. As we will discuss in Sec. II a nonzero 
density of neutrinos always increases fiQ. If the baryon 
chemical potential is small and the chemical potential of 
the neutrinos becomes large enough, it will be energeti- 
cally favored to convert some of the neutrinos into pions 
and electrons. In this way as we will show in Sec. II a 
pion condensate is formed in a dense neutrino gas for 
low baryon chemical potential. This has interesting con- 
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sequences for the phase diagram of QCD at finite elec- 
tron and/or muon lepton number chemical potential. In 
Sec. II and Sec. Ill we will show that pion condensation 
takes place in a large part of the phase diagram of QCD 
for \hb\ m p ~ m ir- Even before the onset of pion con- 
densation, the neutrino gas has interesting implications 
on the behavior of the pions. Because of the nonzero 
isospin chemical potential, the degeneracy in the masses 
of the charged pions will be lifted. Moreover, the lifetime 
of the pions will change as we will see in Sec. IV. 

Let us point out that similar condensation phenom- 
ena can also arise at larger baryon chemical potential 
(/xb > m p ). In dense electrically neutral nuclear matter, 
—\iq can become of the order of the vacuum pion mass. 
Hence condensation of negatively charged pions in dense 
nuclear matter seemed a realistic possibility [23L I24I [25[ . 
However, due to interactions the in-medium pion mass 
in nuclear matter is increased such that pion condensa- 
tion becomes improbable. On the other hand, the mass 
of the negatively charged kaon is decreased due to an 
attractive interaction. A negatively charged kaon con- 
densate turns out to be a more likely possibility in dense 
neutral nuclear matter [2o| even at high neutrino den- 
sities [27]]. In color superconducting matter which is 
formed at even higher baryon densities, the masses of 
the kaonic excitations are smaller than those of the pio- 
nic excitations [28|. It turns out that in electrically neu- 
tral color superconducting matter a neutral kaon conden- 
sate can be formed-creating the so-called CFL- K a phase 

In nature a dense neutrino gas is created during a su- 
pernova explosion. Part of the neutrinos produced in this 
explosion will be trapped in the precursor of a neutron 
star, the so-called proto-neutron star (36|. Inside a proto- 
neutron star the baryon density is very large, which as we 
stated above makes pion condensation unlikely. However, 
one could speculate that in the crust or in the atmosphere 
the neutrino density might be large and baryon density 
small enough for pion condensation to occur. To fully 
settle this issue one should also show that chemical and 
thermal equilibrium can be reached. In this paper, we 
will not give a definite answer to the question whether 
or not a pion condensate is realized somewhere in our 
universe. We will only unambiguously show that a pion 
condensate is formed if the conditions are right, that is 
in a dense neutrino gas at low baryon densities. 
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II. PION CONDENSATION 



potential of the charged leptons is given by 



At very large neutrino density and small baryon chem- 
ical potential it is energetically favored to convert some 
of the neutrinos into pions and leptons. In this section 
we will show that in such a case a pion condensate will be 
formed. We will compute at which neutrino density this 
will happen and discuss the phase diagram as a function 
of electron and muon lepton number chemical potential. 

In order to describe matter at finite density we express 
the chemical potentials of the up quark (u), down quark 
(d), electron (e), electron neutrino (V e ), muon (/x) and 
muon neutrino (i/„) in terms of the chemical potentials 
of conserved quantities. We will assume that chemical 
equilibrium is reached, hence 
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Here hb, HQi Hl c and hl , denote respectively the chem- 
ical potentials for baryon number, electric charge, elec- 
tron lepton number and muon lepton number. Strictly 
speaking, the muon and electron lepton numbers are not 
separately conserved due to neutrino oscillations. We 
will assume that the size of our system is much smaller 
than the neutrino oscillation length. Furthermore, we 
will assume that the time it takes to bring the system 
in chemical equilibrium is much shorter than the typi- 
cal neutrino oscillation time. In this situation it is to 
first approximation correct to introduce separate chemi- 
cal potentials for electron and muon lepton numbers. On 
long time scales only total lepton number is conserved 
and one can no longer talk about separate lepton num- 
ber chemical potentials. This situation is a special case 
of our results in which hl c = hl ■ 

To make the system of particles electrically neutral, 
the charge chemical potential has to be chosen in such 
a way that the charge density uq vanishes. This can be 
achieved by solving the following equation 
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where Q denotes the thermodynamic potential. It has a 
hadronic and leptonic component. We will write 
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where Qh is the contribution of the hadrons, fit that of 
charged leptons and Q, Vi that of the neutrinos. To a good 
approximation the leptons can be described by a free non- 
interacting gas of fcrmions. Hence the thermodynamic 
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where uj p = {p 2 + m 2 ) 1 / 2 and = 1/T denotes the inverse 
temperature. For the neutrinos (which only have left- 
handed chirality) one has 
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The densities of the leptons can be found by taking the 
derivative of the thermodynamical potential with respect 
to the chemical potential. In this way for zero temper- 
ature one finds that the density of the charged leptons 
equals 



m) [(^L, - HQ? - rn\ 
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for \nLi — Hq\ > nii and rii 
of the neutrinos is equal to 



otherwise. The density 
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Now let us set the temperature to zero and take 
Hb small so that there are no baryons in the system. 
This means that the chemical potential of the proton 
(p p = hb + Hq) should be smaller than the proton mass 
m p . Since hq can become m n , this translates into the 
requirement that \hb\ m p ~ m -n- In that case it is 
well known from the analysis using the chiral effective 
Lagrangian that a 7r + condensate will form via a second- 
order transition if hq > ra, [l|, 0] . This is confirmed by 
lattice QCD simulations 

If in this particular situation \hq\ is smaller than m w , 
there are no hadronic states around. Hence the hadronic 
sector is automatically electrically neutral. The question 
then is, can we neutralize the leptonic sector such that 
Hq = TOtt? If so, we have found the onset of pion con- 
densation in electrically neutral matter. Answering this 
question amounts to solving the following equation 



^V;(/iQ = m n ) = 0. 



(8) 



If there arc no neutrinos around this equation has no so- 
lution since hq = automatically introduces a nonzero 
density of electrons and muons. However, if we for ex- 
ample take fj,L B = Hl = Hl we find that HQ = hl guar- 
antees electric neutrality. Pion condensation then sets in 
at hl = m-K- This corresponds to a electron and muon 
neutrino number density of m^/(67r 2 ) s» 5.9 x fm~ 3 . 



A. Phase diagram /il c vs. fiL^ at T = 

Now that we found the onset of pion condensation in 
electrically neutral matter for equal electron and muon 
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FIG. 1: Phase diagram of electrically neutral QCD at T — 0, 
valid for \hb\ ^ rn p — TH-k, as a function of the electron and 
muon lepton number chemical potential. The phase boundary 
denotes a second-order transition. The black dot denotes the 
point where ul^ — 0. 



FIG. 2: Phase diagram of electrically neutral QCD, valid for 
I/is | i$ rn p — m-Tv, for T — (solid, red), 25 (dashed, green), 
and 50 MeV (dotted, blue), as a function of the electron and 
muon lepton number chemical potential. The T — 25 and 
50 MeV results are estimates. 



lepton number chemical potential, let us see what hap- 
pens in the more general case in which they are different. 
By solving Eq. ((SJ) for (1q = we obtain the phase dia- 
gram of QCD as a function of /Lt£ e and for T = and 
small baryon chemical potential. We display the result 
in Fig. HJ 

Let us discuss the transition line to the 7r + -condensed 
phase in somewhat more detail. If m, — < (1l < 
m.n + there are no muons in the system so it should 
be neutralized solely by the electrons. This then gives 
to^ — m e < /Z£ e < + m e which corresponds to the 
vertical lines in the phase diagram of Fig. [TJ Hence the 
length of the vertical line is equal to twice the muon mass. 
For Hl c < ~ m e the exact solution to the neutrality 
constraint gives the following phase boundary 

= rn v + J (fi Lc - m^) 2 ~m 2 e + m£ . (9) 

while for /i£ e > m, + m e we find 

Mz„ = - sJ{pL e -m^) 2 -ml + ml. (10) 

If the muon lepton number chemical potential vanishes, 
a 7r + condensate will form if 

fiL e >m 7r + ^Jml + m 2 e - m 2 ^ = 231 MeV. (11) 

This corresponds to an electron neutrino density of 
mL/(6tt 2 ) = 2.7 x 10- 2 fm~ 3 . 

In some situations like the proto-neutron star it is 
natural to require zero muon lepton number density 
(nt, = rip + n v =0). In such a case we find that 
the transition to the pion condensed phase occurs at 

Mi„ = £m w - ^/m2 +( m 2 - m 2)/£ = 31 MeV, (12) 



here £ = 1/[1 — (1/2) 2 / 3 ]. The electron chemical potential 
is then equal to 

= + ^)%+^. (13) 

To a good approximation we can neglect the electron 
mass in the last equation. Then by expanding Eq. (|12p 
in powers of (m 2 — m 2 )/m 2 we obtain 

(1 m 2 — m 2 \ 
1+ 2^^r^J =164 MeV. (14) 

The solution using the exact expression gives = 
163 MeV. Hence if the muon lepton number density van- 
ishes pion condensation will occur if the electron neutrino 
density is larger than /z 3 /(67r 2 ) = 9.6 x 10~ 3 fm . 

B. Phase diagram \XL e vs. fj,z, for T -C rn.Tr 

Let us now discuss how the fiL c vs. Ml m phase diagram 
of electrically neutral QCD looks like at finite tempera- 
ture. As long as the temperature is much smaller than 
the chiral phase transition temperature T c , the mass of 
the pion is not much different from the vacuum pion mass 
m^. Hence for T <C T c the onset of pion condensation 
still occurs to a good approximation at \iq = m % . 

If the temperature is also much smaller than m v , there 
are to a good approximation no hadronic states in the 
system as long as fj,Q < m^. Hence we can like in the 
previous subsection find the onset of pion condensation 
by solving Eq. (|5J), but now at finite temperature. We 
display the results in Fig. [2] 
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To understand Fig. [2] better, let us see what hap- 
pens if the temperature is increased. In that case the 
muon mass, which is equal to half of the length of the 
vertical line in the phase diagram, becomes less impor- 
tant. Let us for a moment neglect the muon mass com- 
pletely. Then the solution of Eq. ([8|) is a straight line 
from (fJ.L e , = (2m-7r,0) to (0,2771^). At these high 
temperatures where the muon mass can be neglected, the 
approximation of ignoring the hadrons is wrong, so that 
the straight line will never be the real phase boundary. 
Nevertheless, this argument shows that by raising the 
temperature the tendency of the phase boundary is to ap- 
proach the straight line as long as T <C m v in agreement 
with the results displayed in Fig. H Asa result of this 
straightening tendency there are points in the phase dia- 
gram (fiL e < fn-K-, PL > fn n ) where the pion condensate 
does not appear at T — but arises if the temperature 
is increased. 



III. NJL MODEL CALCULATION 

As was shown in the previous section, pion conden- 
sation arises in electrically neutral hadronic matter for 
large neutrino densities. We have obtained the phase di- 
agram of QCD for low temperatures and small baryon 
chemical potential. It is of interest to see what happens 
to the phase diagram at larger baryon chemical potential 
and higher temperatures. 

To investigate the phase diagram at finite T and \ib we 
need to know the hadronic component of the thermody- 
namic potential, f2/j. Unfortunately we can not compute 
ilh in QCD from first principles. We will resort to an 
effective model of QCD, the Nambu-Jona-Lasinio (NJL) 
model (see e.g. Ref. [37| for an extensive review). The 
NJL model qualitatively describes QCD. It for example 
captures the chiral phase transition and can be used to 
study mesons and pion condensation. However, it is im- 
portant to realize that nuclear matter is not described at 
all with this model. The phase diagrams obtained with 
the NJL model will therefore only describe the qualita- 
tive features of the QCD phase diagram. Nevertheless, 
as we argued in the previous section, our conclusion that 
a large part of the phase diagram contains the pion con- 
densate is certainly also correct for QCD. 

The 2-flavor NJL model we will use in this article is 
given by the following Lagrangian density 

C = -0 (i-f^d^ - m + /^yo) tp + (1 - a)Ci + aC 2 , 

A = GpV) 2 + (V^) 2 + (^75^) 2 + (^75^) 2 ], 

C 2 = G[{^) 2 - (VifVO 2 - (^75^) 2 + (^75r^) 2 ]. 

(15) 

Here G is the coupling constant and f denotes the Pauli 
matrices. The mass matrix m is diagonal and contains 
the bare quark masses m u and m d . The matrix p is also 
diagonal and contains the quark chemical potentials p u 
and p d , which are given in Eq. {!]). The interaction terms 



L\ and L 2 arc invariant under SU(2)l x SU(2)r trans- 
formations, as is QCD. The interaction L\ preserves the 
axial U(1)a symmetry, while £2 (the two-flavor 't Hooft 
term) breaks it. In QCD the U(1)a symmetry is broken, 
so a should be nonzero. 

Since the NJL model is an effective model which is only 
valid for low momenta, one has to specify an ultraviolet 
regularization. In this paper we will employ a three- 
dimensional sharp momentum cutoff A. 

We will take the following parameter set which is also 
used in other papers (see e.g. Refs. [H, H^|): G = 

5.04 GeV~ 2 , A = 651 MeV and m = m u = m d = 

5.5 MeV. We will treat a as a free parameter. At 
zero temperature and baryon chemical potential, its 
most likely value lies somewhere between 0.1 and 0.2 
[37l . [40| . With this parameter set one can reproduce 
the vacuum properties of the mesons, m^± « 140 MeV, 
U ~ 92 MeV. 

In order to compute the thermodynamic potential, we 
will resort to the mean-field approximation. We will in- 
troduce the following mean fields 



a u = -4G{uu), 
a d = -4G{dd), 
P = -2G(ipT2i r y5i'), 



(16) 



here a u and ad are proportional to the chiral condensates, 
while the pion condensate is described by p. 

The hadronic part of the thermodynamic potential 
can now be found by expanding the Lagrangian density 
around the mean fields and integrating out the fermions. 
In that way one obtains 



tt h = 

Nc_ 
2tt 2 



|(1 - a)(er 2 + aj) + aa u a d + p 2 
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Here N c = 3 denotes the number of colors and A; are the 
four independent eigenvalues of the mean-field Hamilto- 
nian TL which is given by 



n 



7°7 • p + 70 - Hu 
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where the constituent quark masses are equal to 
M u = m u + (1 - a)a u + aa d , 



M d 



m d 



[1 - a)a d + aa u 



(18) 



(19) 



If M u = M d = M the four eigenvalues can be obtained 
analytically and read 



Aj = yj (u p ± Mq/2) 2 + P 2 ± M 



(20) 



where p, = /is/3 + pg/6 and uo p = \/p 2 + M 2 . The 
values of the mean fields follow by minimizing the ther- 
modynamic potential with respect to these mean fields, 
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FIG. 3: Evolution of the condensates and chemical potentials along the \a,b = axis at zero temperature. Left panel: The 
fiL^ = case. Right panel: The nz^ = case. Displayed are: constituent quark mass (red, solid), pion condensate (green, 
long-dashed), electric charge chemical potential (blue, short-dashed), in-medium ir + mass (magenta, dash-dotted), and muon 
lepton number chemical potential (cyan, dotted). The in-medium pion mass is discussed in Sec. II V Al 



which amounts to solving the following equations 



on _ on _ on 

da u da d dp 



0. 



(21) 



In addition one has to solve the electric neutrality con- 
straint, Eq. and in the case where we require zero 
muon lepton number density, also the constraint til =0. 
If a = 1/2, the effective potential is a function of a u -{-as- 
Hence in that case one always finds a u = a d as a solution. 
For other values of a, a u and as, can be different. 

We have solved Eq. (|2"Tj) and the neutrality constraints 
numerically and will discuss the results in the following 
subsections. 



A. Condensates at T = and /is = 

In Fig. [3] we show the results achieved at zero temper- 
ature and baryon chemical potential. In this limit the 
condensation is governed by the chiral physics and the 
NJL model, augmented with the gas of free leptons to 
account for charge neutrality, is therefore expected to be 
most reliable. 

First of all, note that the NJL calculation naturally 
reproduces the previous model-independent evaluation of 
the onset of pion condensation. To demonstrate that 
even the quantitative results inside the pion-condensed 
phase are to a large extent model- independent, one can 
employ chiral perturbation theory [l|, Q. The leading- 
order (Euclidean) chiral Lagrangian reads 

£ xPT = f -f [TriDprfDpU) - 2m%ReTrU] , (22) 



where U is a unitary 2x2 matrix field and the covariant 
derivative, D^U = d^U—i^U +Ur fi , incorporates its cou- 
pling to external left- and right-handed vector fields, 
and . In presence of electric charge or isospin chemical 



The ground-state expectation value of the order pa- 
rameter can be without lack of generality cast as 



U 



cos 8 sin 8 
- sin 8 cos 8 



(23) 



If 8 = we recover the vacuum chiral condensate, while 
8 = it/2 would be a purely pion condensate. In general 
this rotation angle is determined by minimization of the 
static part of the Lagrangian and for \/J,q\ > is given 

by H 



cos e = 



(24) 



This is analogous to the gap equation in the NJL model 
and provides a model-independent relation among the 
chiral and pion condensates and the electric charge chem- 
ical potential. The electric charge density of the pion 
condensate is in turn obtained as [ij 



dC 



X PT 



dfiQ 



(25) 



potential, they arc 



We used these expressions supplemented with the lep- 
ton sector to check the results of Fig. [3] Since chi- 
ral perturbation theory allows a straightforward evalu- 
ation of the expectation values of fcrmion bilinear oper- 
ators rather than the NJL order parameters cr Uj d, p, see 
Eq. (|16[) . we use the vacuum NJL value of a u = as to 
normalize the condensates obtained with chiral pertur- 
bation theory. The results agree with NJL model up to 
about a 2% accuracy even at the highest values of hl c 
considered. Therefore, we do not plot the correspond- 
ing lines in Fig. [31 they would hardly be distinguishable. 
(See fl3j for a similar comparison.) 
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FIG. 4: Phase diagram of the electrically neutral NJL model 
with a = at various values of temperature (indicated by the 
numbers in boldface). The black dots denote the points at 
the different phase boundaries where the muon lepton number 
density (nz,„) vanishes. 



B. Phase diagram \Xh e vs. \it, at finite temperature 

In Fig. [4] wc display the phase diagram of the NJL 
model as a function of fj,L e and [1l f° r different tem- 
peratures at \ib = and a = 0. We found that up to 
about T = 150 MeV this phase diagram does not change 
much by modifying a. Comparing this phase diagram 
to Fig. [2] we find very good agreement with the model- 
independent estimates up to T = 50 MeV. Hence we can 
say that Fig. [2] most likely represents the phase diagram 
of electrically neutral QCD for small baryon density as a 
function of lepton number chemical potentials. 

For higher temperatures the NJL model calculation 
serves as an illustration of how the phase diagram of QCD 
would qualitatively look like. We find that the onset of 
the pion condensed phase moves to higher values of the 
lepton number chemical potentials. This is caused by 
the appearance of charged hadronic states in the system 
which affects the neutrality. Furthermore due to the in- 
crease of the pion mass at finite temperature, the onset 
of pion condensation moves to larger values of /iq. 



C. Phase diagram (ib vs. (il at zero temperature 

In Ref. [4l| the phase diagram of neutral matter with 
neutrinos present was discussed for very large baryon 
chemical potential where color superconducting phases 
arise. See also Refs. [H, [43|. Here we complement this 
phase diagram with the results obtained at low baryon 
chemical potential. We display the results in Fig. [5] for 
the two cases regarding the muon-lepton content dis- 
cussed above, and for two different values of the a pa- 



rameter. 

The structure of the phase diagram is simple. The 
straight line, which determines the onset of pion con- 
densation at low baryon chemical potential, is given by 
the model independent argument presented in Sec. II. 
The corner of the pion condensed phase, where this line 
breaks, is marked by the appearance of up quarks in 
the system. In the NJL model at the mean-field ap- 
proximation, quarks behave as noninteracting quasipar- 
ticles with effective energy gap (mass) determined by 
Af«,eff = M u - \[jt B - M d , e s = M d - + i/iQ. 

When this drops below zero, a Fermi sea of quarks is 
formed. This happens along the boundary of the pion 



condensed phase, when ^ = M u 



Using the con- 



stituent quark mass in the vacuum, M u d = 325 MeV for 
our parameter set, we find = 232 MeV. 

The line of second order phase transition eventually 
ends up in a critical point from which on the transition 
becomes first order. To demonstrate this, we plot in Fig. 
[6] the values of the condensates, the charge chemical po- 
tentials, and the in-medium pion mass (obtained as the 
pole of the pion propagator [2(|) along the section of the 
phase diagram with fj,L e = 400 MeV for the =0 case 
and a = 0. We observe that as a decreases, the critical 
point moves down and the pion-condensed phase shrinks 
to somewhat smaller [ib- This is because small a tends 
to split the constituent masses of up and down quarks 
and thus disfavor the pairing mechanism which underlies 
the pion condensation. 



IV. PION PROPERTIES IN A NEUTRINO GAS 

As we saw in the previous sections, the lepton medium 
induces an isospin chemical potential. This chemical po- 
tential and the finite density of neutrinos will modify the 
behavior of the pions. Let us therefore have a closer 
look at the spectral properties of the pions in the lepton 
medium, in particular their masses and decay rates. 



A. Masses 

As long as the isospin is not spontaneously broken (i.e., 
there is no pion condensate) , the in- medium pion masses 
are simply given by M n ± — =p /iq [l|, Q • The pion 
condensation sets in where one of the masses drops to 
zero. 

At low jU£ e , only electrons are light enough to be ex- 
cited. In order to preserve electric neutrality, we have to 
make sure that the electron chemical potential is zero, 
that is, /Iq — /U£ e . The pion masses are thus simply 
equal to M n ± = m v =F hl c - When fj,L s exceeds the muon 
mass, neutrality becomes a nontrivial issue and the pion 
masses depend on the exact way it is imposed. We will 
describe in detail two special cases: Zero muon lepton 
number chemical potential, and zero muon lepton num- 
ber density, both of which were already discussed in Sec. 
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FIG. 6: Evolution of the condensates, chemical potentials 
and the in- medium pion mass as a function of fiB for (j,L e — 
400 MeV, a = and T = in the n L =0 case. 



II. 

At fjLi, = we have /i„± = ±/iq so that for hl c > 
rn M + m e , antimuons will appear in the system. In this 
case electric neutrality requires that the electrons and 
antimuons have the same Fermi momentum, which leads 
to 



fJ-Q 



A. 



(26) 



Setting fiQ = m n recovers the transition point to the 
pion condensed phase, Eq. (11). 

If we instead demand zero muon lepton number den- 
sity, the /i + chemical potential modifies to pq — (J-l^ and 
fiL is determined self-consistently from the muon lep- 
ton number neutrality condition; a negative contribution 
from antimuons has to be compensated by a finite den- 
sity of muon neutrinos. Solving the set of two neutrality 
equations, we obtain the results shown in Fig. [3] by the 



dash-dotted lines. The mass of 7r is given simply by a 
reflection of the curves with respect to the M = line. 

B. Decay rates 

The charged pions in the vacuum decay predominantly 



in the 7t + 



channels. These consti- 



tute about 99.99% of the total decay rate [44(. In the 
vacuum and at rest, the decay rate is, at the leading or- 
der in weak interactions, given by the textbook formula 



1 

4- 



V ud G Ffl m 



1 



(27) 



Here V denotes the Cabibbo-Kobayashi-Maskawa ma- 
trix and Gf the Fermi constant. The presence of the 
medium breaks explicitly Lorentz invariance, and we will 
therefore calculate the decay rates as a function of mo- 
mentum. The derivation entails two ingredients which 
can be considered separately: The invariant decay am- 
plitude and the kinematics. 

To determine the invariant amplitude, we need to know 
the coupling of the pion to the charged weak current. In 
vacuum, this is equal to ^gVudj-nV^ > where p M is the pion 
four-momentum and g the weak coupling constant. To 
determine how this is modified in the medium, we use 
chiral perturbation theory [45| . Introducing in the lowest 
order pion effective Lagrangian (|2"2"|) the external vector 
electromagnetic as well as the left-handed charged weak 
current, we find that one simply has to make the replace- 
ment — ► p M = (po + pQ,p). In this simple estimate 
of the decay rate, we neglect truly quantum corrections 
to the pion decay constant which give it a weak medium 
dependence (4d| . 

An explicit calculation shows that this change of the 
pion coupling to the weak current is precisely what is 
needed to make the total amplitude for the pion de- 
cay completely independent of the chemical potentials; 
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FIG. 7: Decay rate of n + in units of the vacuum rate at rest 
as a function of the pion velocity for various values of electron 
lepton number chemical potential (indicated by the numbers 
in boldface; the same values are implied in the lower plot). 
The thin dotted line indicates the vacuum decay rate, i.e., is 
trivially given by the 7-factor. Upper panel: The = 
case. Lower panel: The til u = case. 



the extra term in the pion-weak-current coupling can- 
cels with a similar term coming from the muon and neu- 
trino wave functions. As a consequence the kinematics 
becomes trivial: The integration over phase space for the 
final state can be transformed into the center-of-mass 
frame (CMS). The only effect of the medium is a restric- 
tion of the phase space due to the fact that some of the 
momentum states are occupied by particles in the Fermi 
sea. The final formula for the decay rate as a function of 
momentum p then reads 




1 UJ* 

Ip 47r' 



(28) 



where 7^ denotes the Lorcntz factor and to* is the allowed 
solid angle for the products of the decay, measured in 
CMS. 

First of all, observe the decay rate of ir~ is trivially the 
same as in vacuum, because there are neither muons nor 
muon antincutrinos in the system. In the following we 
will therefore concentrate on the decay of ir + . In CMS 
the magnitude of the momentum of the antimuon and 



muon neutrino from the 7r + decay is 



fc 



2m x 



(29) 



the decay is isotropic. (Note that since both the electric 
charge and the lepton number are preserved in the decay, 
the presence of chemical potentials docs not affect the 
energy and momentum conservation.) Therefore, when 
the Fermi momentum of both antimuons and neutrinos 
is smaller than fc , the decay rate will be r , otherwise it 
will be zero, at least for the considered process. 

The energies of the antimuon and the neutrino in the 
rest frame of the medium are given by a simple Lorentz 
boost, 

£/i =7p(eo + /%fcocos6>*), e v =7^0(1 - /%cos0*), 

(30) 

where j3g denotes the velocity. Furthermore, eo = (m^ + 
mj l )/(2m v ) is the antimuon energy and 9* the angle of 
the antimuon with respect to the pion momentum, both 
measured in CMS. This allows a straightforward deter- 
mination of the available phase space for the decay. For 
instance, in the case /i£ = 0, only phase space blocking 
by antimuons occurs (there arc no muon neutrinos) and 
one may explicitly express to* as 



47T 



1 



f3pk V 7^ 



e 



(31) 



where the value of the function S(x, a, b) with a < b is 
equal to a if x < a, to b if x > b, and to x otherwise, that 
is, S(x, a, b) = min[max(a;, a), b\. 

The numerical results for the decay rate as a function of 
the pion velocity are shown in Fig. [Jj The interpretation 
of the results is simple. If the chemical potentials are low 
enough such that the decay is not blocked in CMS, then 
boosting to finite momentum may bring the backward 
emitted particles into the Fermi sea and thus suppress the 
decay rate beyond the simple 7-factor from time dilation. 
On the other hand, if the decay is Pauli-blocked in CMS, 
boosting to finite momentum may liberate the forward- 
emitted particles so that the decay becomes possible. In 
the riL =0 case parts of the phase space are blocked by 
both antimuons and muon neutrinos. 



V. CONCLUSIONS 

In this article we have shown that a positively charged 
pion condensate arises in electrically neutral matter at 
high neutrino and small baryon densities. We found that 
at zero temperature, zero baryon chemical potential, and 
zero muon lepton number density, the onset of pion con- 
densation lies at an electron neutrino number density of 
9.6 x 10~ 3 fixT 3 . 

For zero temperature we have obtained the phase 
diagram of electrically neutral QCD as a function of 
/tL e and /xl , valid for small baryon chemical potential, 



9 



MbI ^ m p ~ m 7r- We have estimated this phase diagram 
at higher temperatures. Comparison to model calcula- 
tions with obtained using the NJL model show that up 
to about T = 50 MeV these estimates are very good. 

Using the NJL model we also studied the behavior of 
the phase diagram as a function of baryon chemical po- 
tential and lcpton number chemical potential. We found 
that the pion condensed phase arises up to |/is/3| ~ 
A I — ^TOtt, where M is the constituent quark mass in 
the vacuum. 

Hence the main conclusion of this work is that the pion 
condensed phase makes up a large part of the phase dia- 
gram of electrically neutral QCD at finite lepton number 



chemical potential and small baryon chemical potential. 
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